Comments on "Analytic and Numerical Study of Preheating Dynamics" 
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In a recent paper Boyanovsky et al ( |iep-ph/9608205 ) studied preheating for a particular theory of 
the inflaton field. Although they results essentially confirm those of the original papers on reheating 
after inflation, they claimed that most of the earlier works on preheating did not take into account 
backreaction of created particles, misused the Mathieu equation, and were inconclusive about the 
symmetry restoration after preheating. We explain why we cannot agree with these statements. 
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Recently there was an important progress in the the- 
ory of reheating of the universe after inflation. As it 
was shown in in many versions of chaotic infla- 

tion the inflaton field very rapidly decays into its own 
quanta and other bosons due to broad parametric reso- 
nance. To distinguish this stage of explosive reheating 
from the stage of particle decay and thermalization, we 
called it preheating. Investigation of this process should 
be performed both in the broad and in the narrow res- 
onance regime, with an account taken of backreaction 
of created particles and of expansion of the universe [Q . 
Bosons produced at that stage are far away from thermal 
equilibrium and have enormously large occupation num- 
bers. Preheating leads to many interesting effects. For 
example, specific nonthermal phase transitions may oc- 
cur soon after preheating, which are capable of restoring 
symmetry even in the theories with symmetry breaking 
on the scale ~ 10^^ GeV [§,1). Strong deviation from 
thermal equilibrium and the possibility of production of 
superheavy particles by oscillations of a relatively light 
inflaton field may resurrect the theory of GUT baryoge- 
nesis ||^ and may considerably change the way baryons 
are produced in the Affleck-Dine scenario jsj. 

The theory of preheating is rather complicated and 
strongly model-depended. Certain aspects of this theory 
were investigated by many authors, and a lot of interest- 
ing information has been accumulated [p^ . However, 
recently Boyanovsky, de Vega, Holman, and Salgado have 
written a paper criticizing almost all other authors study- 
ing preheating [ p7[ . They write, in particular: "There 
have been a number of papers (see refs. [0j|J|,§|) dedi- 
cated to the analysis of the preheating process... In most 
treatments this approximation neglects all back reaction 



effects. We will argue that these approaches do not agree 
with our analytical and numerical results and miss im- 
portant physics... In fact, whereas the Mathieu equation 
has infinitely many forbidden bands, the exact equation 
has only one forbidden band. Such an approximation, 
based on the chart of the unstable bands of the Mathieu 
equation has been used in most etc. treatments of 

preheating. These results must therefore be regarded as 
highly suspect... We also discuss why the phenomenon of 
symmetry restoration at preheating, discussed by various 
authors g- H is not seen to occur." 

We believe that these statements, made in in bold 
face and in italic, are incorrect. They contradict the con- 
tents and misinterpret the papers they cited. 

1. Back reaction effects 
Creation of particles leads to several effects which can 
change the dynamics of the system. Already in our first 
paper [Q we did take into account backreaction of cre- 
ated particles, which modifies equations of motion for the 
interacting fields. We also took into account the contri- 
bution of the created particles to the energy density of the 
universe, which controls the speed of the universe expan- 
sion. Khlebnikov and Tkachev studied the feedback 
effects of rescattering of created particles. A section of 
the review article by Kofman in Q| is devoted to the dis- 
cussion of all of these feedback effects of created particles. 
Ironically, only part of the feedback effects (modification 
of equation of motion without an account taken of ex- 
pansion of the universe and rescattering) was considered 
in|l7|. 



2. The Mathieu equation 
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The main subject of investigation of ref. Q was the the- 
ory of a massive inflaton field interacting with another 
scalar field x- Parametric resonance in this theory is 
described explicitly by the Mathieu equation, and com- 
ments of ref. simply do not apply to this theory. 

On the other hand, the main subject of investigation of 
[ prt was the theory m^(/)^/2 -I- Xcp'^/A. There are two main 
regimes there: (f> <C m/\/A and cj) ^ m/V\. The authors 
of do not make any distinction between these two 
regimes because they neglect expansion of the universe. 
However, from [Q it follows that in expanding universe 
there is no parametric resonance at all for cj) m/\/\. 
Therefore all results of related to parametric reso- 
nance in this regime do not give a correct description of 
reheating in the theory rr?<^'^ jl -I- A0'*/4. 

As for the regime 3> m/^/X, the main statement of 
ref. is that the parametric resonance in this model 
is described by Lame equation rather than by Mathieu 
equation. However, this fact was already emphasized ear- 
lier in j^]. To study parametric resonance in this theory, 
one can approximate Lame equation by Mathieu equa- 
tion. The resonance which was found in 0| for this par- 
ticular theory was only in one band, so the statement of 
ref. ijl^ that we used charts of many bands is misplaced. 

Approximate treatment of Lame equations in terms of 
Mathieu equation was made in 1^ only as a first illus- 
trative part of a complete investigation. It was found 
that fluctuations approximated by Mathieu equation 
the unstable modes grow as e^ '*^^'^*, where /i ~ 0.07. 
Direct investigation of the Lame equation reveal similar 
result but with /i ~ 0.04, see the review by Kofman in [Q. 
This was the value which was used in our subsequent pa- 
per Q . Khlebnikov and Tkachev [Q made a more com- 
plete numerical investigation of preheating in the theory 
A(/)''/4 than the investigation performed in jl^, because 
they took into account rescattering of produced particles 
and expansion of the universe. They obtained the same 
value of /i All of these results were ignored by the 
authors of ref. [ pTf . 

A more general theory studied in [Q| contains both the 
terms r7i^0^/2 -|- A0^/4 and g'^4>^x^- In the hmit (j) < 
m/y/\ decay due to parametric resonance can go only 
to the field Xi and there are many instability bands in 
this regime. As for the regime 3> to/a/A, S> A, the 
number of instability bands for the field x is infinitely 
large for all but very special relations between and A. 

It is instructive to consider the maximum value of 
growth index /j, in the unstable bands as a function of 
g^/A. It turns out that the model considered in [ p7| in a 
certain sense is a degenerate case. The function fl(g^ / X) 
has a minimum at g^ / X — 3. In this case equation de- 
scribing creation of x-particles is similar to the equation 
for 0-particles in the theory A(/)*/4. In this particular 
case there is a single instability band. However, in the 
regime 3> A there are many instability bands, and the 
results for are converging to those obtained from the 
Mathieu equation. 



3. Symmetry restoration from preheating 
Finally, let us consider the theory of nonthermal phase 
transitions ||^,^ . These phase transitions have been found 
in several versions of chaotic inflation theory when the 
inflaton field rolls down from very large values (j) Mp. 
Then the oscillations of this field produces enormously 
large fiuctuations {(jP) and (x^), which are capable of 
restoring symmetry in the theory. 

The authors of ref. (l^ claim that they did not find any 
symmetry restoration from preheating. However, they 
studied a regime which is completely different from that 
of 1^,^ . Instead of analysing fluctuations produced by the 
field rolling from (j> ~ Mp, they studied the theory where 
the field (f> falls to the minimum of the effective potential 
from the vicinity of = 0. Symmetry restoration in 
this case would imply return of the field (j) exactly to the 
point (p — 0, which is hardly possible because it looses 
some part of its energy for particle production. Thus it is 
not surprising that Boyanovsky et al did not find in their 
model the effect of symmetry restoration discovered in 
&|6l in a different context. 
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